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Abstract 

This paper presents an approximation algorithm for finding a shortest path between two points s and t in 
a weighted planar subdivision V. Each face / of V is associated with a weight Wf, and the cost of travel 
along a line segment on / is Wf multiplied by the Euclidean norm of that line segment. The cost of a path 
which traverses across several faces of the subdivision is the sum of the costs of travel along each face. Our 
algorithm progreeses the discretized shortest path wavefront from source s, and takes polynomial time in 
finding an e-approximate shortest path. 

1 Introduction 

Let V be a planar subdivision, specified by n vertices and 0(n) faces. We assume, w.l.o.g., that the faces are 
triangles. Each face / has a weight, aj £ Z + , associated with it. This weight specifies the cost per unit distance 
of traversing the face /. For simplicity, we assume that the cost of traversing a unit distance along an edge e, 
say a e , which bounds faces fi and f'2 is mm{a.f 1 ,a.f 2 }. The cost of a straight line path between two points p 
and q on face / is a/||p#||. And the weighted length of a path passing through a number of faces is the sum of 
the weighted lengths of its sub-paths through each face. Let s and t, respectively, be the source and destination 
points on V. The weighted shortest path problem is to find a path having the minimum weighted length from s to 
t. The path we need to find is geodesic in the sense that it is locally optimal and cannot, therefore, be shortened 
by slight perturbations. An optimal path is a geodesic path that is globally optimal. A corresponding problem 
can also be defined when the domain, V is the surface of a polyhedron. The applications of this problem are rich 
and varied, ranging from motion planning, to graphics, injection molding and computer-assisted surgery (3j|7l. 

While there has been previous work on this problem (see Table [[J, a polynomial solution which tackled the 
problem in a systematic way was designed by Mitchell and Papadimitriou Q. Their algorithm uses Snell's Law 
and a continuous Dijkstra method to give an optimal-path map for any given source point s. Because of numerical 
precision issues, the algorithm [7] finds an e-approximation to an optimal path. In this paper too, we consider 
e-approximate paths i.e., the weighted approximate shortest path from s to t is of length at most (1 + e) times the 
weighted distance between s to t of an optimal path. 

Given the higher order polynomial solution for this problem, approximation algorithms are provided by 
Agarwal et al. (T), Cheng S.-W. et al. H, Aleksandrov et al. O, and Reif et al. |EJ|9l whose time complexities 
rely on geometric parameters that are outside the logarithm. However, the time complexities of these algorithms 
are pseudo-polynomial (running time is polynomial in the numeric value of the input i.e., exponential in the 
length of the number of bits in the input). The algorithm by (H has a smaller dependence on n as compared 
with but a higher dependece on e and p (defined in Table []]). Previous algorithms by Lanthier et al. ||5l and 
Aleksandrov et al. Q took a discretization approach by introducing Steiner points on each edge. A discrete 
graph is constructed by adding edges connecting Steiner points in the same triangular region and an optimal path 
is computed in the resulting discrete graph using Dijkstra's algorithm. To avoid incurring high time complexity, 
they have used a subgraph of the complete graph in each triangular region. 

In this paper, we design an efficient polynomial time approximation algorithm for finding e-approximate 
paths in weighted planar subdivisions. Previous schemes for this problem either discretize the space using Steiner 
vertices O or partition the space around a point into cones [6] and process all the discrete objects created. Our 
approach discretizes the wavefront in the continuous-Dijkstra's approach by initiating a set of rays from s, and 
the rays are generated uniformly around a source point. Since the approach of maintaining the wavefront exactly 
is complicated, we do not use it in this paper. Rather, in our approach, we utilize only a subset of rays to guide 
the progress of the wavefront as it expands from the source. When two rays together yield an e-approximate 
shortest distance from s to a vertex v, we initiate another discretized wavefront from v. Thus our propagation 
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Table 1: Comparison of Weighted Shortest Path Algorithms 



processes discrete wavefronts from every vertex and can also be used to solve the all-pairs shortest path problem 
in the domain. This method is also applicable in finding geodesic shortest path between two given points such 
that the shortest path lies on the surface of a 2-manifold whose faces are associated with positive weights. Our 
solution uses several of the Lemmas and Theorems proved in Mitchell and Papadimitriou 0, and extends their 
solution to obtain a better bound. 

The current best known algorithm Q to solve this problem requires 0(K-^ lg (^p) lg j lg \) time, where 

K captures geometric parameters and the weights of the faces. This paper describes an algorithm with 0(n 6 lg(p)+ 
n 4 lg(^)) time complexity, where e' = n ( n -2+ LW } » n is the number of vertices, w and W are the minimum and 
maximum weights among all the weights associated with the faces of V, and L is the maximum length among all 
edge lengths of V. The time complexity of our algorithm is better than the known polynomial time algorithms, 
and it does not contain terms that rely on geometric parameters outside the logarithm. Table Q] compares time 
complexity of our result with the previous work. We assume that all parameters of the problem are specified by 
integers. All the weights are assumed to be greater than or equal to one. Note, that as established in [?!], the num- 
ber of points required to define the interaction of the shortest path map with the domain is 0(n 4 ). Consequently, 
designing a polynomial time algorithm of complexity better than 0(n 4 ) appears challenging. The algorithm Q 
is quadratic w.r.t. these events. Our algorithm takes first steps to provide a sub-quadratic (as a function of the 
number of events) solution. The presented technique is also applicable for the case of finding an e-approximate 
shortest path over the surface of a polyhedron with weighted faces. 

Section |2] provides definitions and properties required for the algorithm. Section [3] lists the details of the 
algorithm with pseudo-code. The conclusions are given in Section [5] Detailed proofs and more figures are 
mentioned in Appendix [6] 

2 Definitions and Properties 

Consider a planar subdivision, V, with triangular faces. Two faces f,f'EV are defined to be edge-adjacent if 
they share a common edge e. 

Definition 2.1 A sequence of edge-adjacent faces T = /i, /2, fk+i ' s called a face-sequence if for every i 
in [l,k], fi is edge-adjacent to /j+i, with edge e, common to both and fi+\. And, the sequence £{T) = 
ei, &2-> ••■) e k i s termed the edge sequence of T. The root of the face sequence T is the vertex v\, when f\ is 
specified by the triangle {vi,V2,v%) where e\ = (i>2, ^3). 

Definition 2.2 Given a face f and one of its edges e, a locally f-free path to x € e is defined to be a geodesic 
path P from s to x such that the path does not pass through B(x, 5) Pi /, V<$ > 0, where B(x, 5) is the open ball 
of radius 5 centered at x. 

Intuitively, a locally /-free path to a point x € e strikes x from the exterior of face / and is locally optimal. 



Background results 

The following lemmas are from Q : 

Lemma 2.1 Let P be a geodesic path on V. Consider any edge e = / n /', where f and f are two adjacent 
faces ofV. Let a e = min{af,afi}. If P intersects the interior of e at point q, known as crossing point, then p 
satisfies Snell's Law satisfying a f sin 6 = a/' sin^', where and 9' are the angles made by P w.r.t. the normal 



at q in the faces f and f, respectively. If ' a e = cef and P shares a segment (a, b) with e, then P satisfies local 
criticality condition striking efrom the side of f at a with an angle of incidence 9 c (f, /') = arcsin(a///a/) and 
exits at b into face f at an angle —9 c (f, /'). The segment (a, b) is known as critical segment. The closer of the 
two points {a, b} to s is known as critical point of entry of path P from face f and the other is known as critical 
point of exit of path P into face f. 

Lemma 2.2 The general form of weighted geodesic path is a simple (not self-intersecting) piecewise-linear path 
which goes through an alternating sequence of vertices and edge sequences such that the path along any edge 
sequence obeys Snell's Law at each crossing point and obeys local criticality condition at each critical segment. 
Further, between any critical point of exit and the next critical point of entry along the path, there must be a 
vertex. 

Lemma 2.3 Let p be a shortest locally f-freepath. Let p' be a sub-path of p such that p' goes through no 
vertices or critical points. Then, p' can cross an edge e at most 0(n) times. Thus, in particular, the last edge 
sequence of p contains each edge 0(n) times, implying a bound of 0(n 2 ) on its length. 

Lemma 2.4 There are at most 0(n) critical points of entry on any given edge e. 
Our Framework 

The algorithm uses rays to span the space V. We denote a ray by the ordered pair r = (u, v), where u is the 
starting point of the ray and v is a second point that defines the ray. For a ray r = (p, w), the point p is termed 
as origin(r) and r is said to have originated from p. Let 1Z(p) denote the set of rays that are originated from a 
point p. 

Our algorithm simulates a discretized version of a continuous wavefront approach. We discretize the wave- 
front using a set of rays 1Z(s) that start from the source and fan out in a cone around the source. When a shortest 
path to a vertex is determined, rays are also initialized from that vertex, and create the set 1Z(v). Each ray when it 
strikes an edge of a face of V, either refracts at the edge or critically reflects. As stated in Lemma l2Tl the portion 
of the polygonal edge that critically reflects the rays (approximating the wavefront) has a critical segment with a 
critical point of entry. Consider a ray r c 6 H{v) that traverses over face / and strikes an edge e = / D /' at a 
critical angle 0(f,f). Note that the set of rays that strike at an angle greater than that angle are not part of any 
shortest paths. A shortest path that uses r c may exit edge e back onto face / emanating from e at the same critial 
angle 8(f, /'). We term the point, p(e c ), at which r c strikes e as the critical point of entry. We also term the 
entire segment of e, from p to an endpoint of e which contains points from where the shortest path comprising r c 
may reflect back, as a critical edge segment, from The sets comprising the vertices of V, critical points of entries, 
and critical edge segments are denoted by V, C, and K respectively. 

Importantly, note that not all the rays 1Z(s) denned above are generated during the algorithm; when a ray 
r £ is actually generated or processed, then the ray r is termed traced. 

To achieve an e-approximation, apart from the rays that have originated from vertices and by critical reflec- 
tions from critical edge segments, the algorithm also introduces additional rays to fill in the gap created by the 
critical reflection. These edges atre generated uniformly spaced in a cone with the critical point of entry as the 
apex of the cone. For a source v, consider the rays from v that strike an edge e = (it, w) from face / at a point p 
and let r be the last ray that refracts onto face /'. Further let r' be the first ray that critically reflects and lies on 
(p, w). Let 6 be the angle that r makes w.r.t normal at the point of strike on edge e and 9' the angle of refraction. 
The discretization results in no refracted rays that make an angle w.r.t. the normal at e that lies between 9' and 90 
degrees. To account for these, rays are generated from p that lie in the cone defined by apex p and extremal rays 
r' and (p, v). The cone lies on the face /'. These additinal rays are termed Steiner rays generated from critical 
points. See Fig. [2] 

The organization of Rays 

The set of rays originating from a vertex v and steiner rays originating from successive critical points that arise 
when a ray from v strikes an edge at a critical angle, can be arranged in a tree structure Tr(v), known as tree of 
rays. The root node of Tr(v) corresponds to the vertex v. All the other nodes of the tree correspond to critical 
points of entry as follows: an internal node w in the tree corresponds to a critical point of entry generated by a ray 
from parent(w) striking an edge e G V at the critical angle. Steiner rays are generated from this critical point 
as described above. Thus with each node u of the tree is associated a set of rays, denoted by 1Z(u), generated 
from u and arranged in angle-sorted order. The tree Tr(v) thus comprises the set of rays (traced or otherwise) 



e 4 

Figure 1 : Tracing rays from a vertex 
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Figure 2: Tracing critically reflected and Steiner rays 



originating from v and internal nodes corresponding to the critical points of entry. Note that all the rays in Tr{v) 
are not necessarily traced. 

A path in the tree Tr{v) corresponds to a sequence of points po(= v),p\,p2, ■ ■ ■ ,p m in V, and a sequence 
of rays tq, ri, . . . , r m where each point p, is a critical point and the origin of a set of rays, and n is a ray that 
originates from pi and strikes an edge of V at p^+i. For a point pj in this sequence, the sequence of points 
{POjPij ■ ■ ■ >Pj-i} are termed as ancestors ofpj (or rj) and points in the set {pj + i,pj+2, . . . ,p m } are termed as 
descendants ofpj (or rj). 

Suppose v G V |J C is either a vertex or a critical point of entry. Consider the rays TZ(v) that originate from 
v and are ordered by the angle they make w.r.t a coordinate axis with origin at v. We define two rays n, Ti that 
originate from v as successive/adjacent whenever r\,ri are either adjacent in the ordered set TZ(v) or are the first 
and last rays in lZ(v). See Figs. Q]and[2l 

Let k G K, be a critical edge segment. Since the shortest path can be reflected back from from any point on 
the critical edge segment, the algorithm will generate rays that start from various points on k. These rays reflect 
back into the face bordering k making an angle with the critical segment equal to the critical angle, as specified 
by Snell's Law. Let 1Z(k) be the ordered set consisting of rays that originate from k. The ordering of rays in 
TZ(k) is achieved by the sorted order of the Euclidean distance of their origins from the critical point of entry, cs, 
corresponding to k. We define two rays n, originating from k to be successive/adjacent whenever origin(r\) 
and origin(r2) are adjacent in the ordering along k. See Fig. [3] 

Consider a ray r that has originated from either w G Tr{v) or from a point p for p G k and k G K.. 
Let £ (r) be the sequence of edges intersected by r as the ray propagates via a sequence of refractions. Let 
e = (vi,Vj) be an edge in £(r) and let £(r,e) be the initial subsequence of edges in £(r) that ends at e. 
We let R(e,v) be the set of rays in Tr(v) that intersect e and let R(e,v,r) = {r'\r' G Tr(v) and either 
£(r', e) is a subsequence of £ (r, e)} or the reverse is true. Consider the set of intersection points of rays on 
e, i.e. I(e,v,r) = {p(e,v)\p(e,v) = r n e,Vr G R(e,v,r)}. Let d(p,q) be the Euclidean distance between 
p and q. Define r' to be the ray such that r' = wgmm T€R ( e v r \{d(p(r), Vi)} i.e., r' G R(e,v,r) is the ray 




Figure 3: Critically reflected rays from a critical segment 



which intersects e at the point pi such that among all points in J(e, v,r), pi is closest to Vi. Similarly, let 
r 'sib = ar § mm re-R(e,t),r){^(p( r )i v j)}- The rays r' and r' sib are termed as siblings and are computed by the 
relation r' sib = S(£(r, e),v). These rays will be traced. See Figs. [T]and[3] All other rays in R(e, v, r) are known 
as intermediate rays between sibling rays r' and r' sib . Note that two sibling rays can be successive rays. 

Let r be a Steiner ray that has originated from a critical point in a tree Tr{v), say at node w € Tr{v) and 
r' be the sibling of ray r, i.e r' = S(£(r, e),v), that originated at node w' £ 7r(v). Let v' be the critical point 
of entry at the least common ancestor of w and w in Tr(v) and let the ordered set S(r) = {ro, r\, . . . , ru = r} 
comprise the traced rays that originate from the sequence of points v' ,v\, . . . ,Vk = w, respectively, in Tr(v). 
The tree path P from v' to w, corresponding to the sequence of critical points and sequence of rays in S(r), is 
termed as the critical ancestor path of r. Note that the critical ancestor path of a ray r is always denned with 
respect to its sibling. 

We define two rays r\ , T2 to be parallel if (i) r\ and T2 traverse the same sequence of faces F and (ii) V/ € F, 
the line segments / Pi r± and / n are parallel. 

With every ray we associate a classification, sourcetype, from the set {vertex, critseg, critptentry}. Here, 
vertex, critseg, critptentry flags indicate that the ray r has originated from a vertex source, reflected from a critical 
segment, or originated as a Steiner ray from a critical point of entry, respectively. 

We will use d w {u, v) to be the length of the shortest weighted geodesic path from u to v. 

Algorithm Outline 

We outline the algorithm, which uses essentially a discretized wavefront approach. The algorithm starts with 
a set of rays R(s), that have the source s as the origin and are uniformly distributed in the ball around s. As 
(approximate) shortest paths are discovered to vertices in V, these vertices are used as origins to initiate rays, that 
are used to determine the shortest path map w.r.t. these vertices. 

Let it be a vertex from which rays are originated. Let F be the set of faces incident to u. For each face f € F 
and for each e G /, the pairs of sibling rays in i?(e, u) are traced further. At any point in the algorithm, there 
will be a number of such sibling pair of rays. The events corresponding to the intersection of these rays with the 
edges of V are ordered using event queues. 

In general, consider the set of rays, R(e±, w). Let e\ = / n /'. Suppose that one of these rays, r, strikes edge 
e\ by traversing along the face /' and is required to be refracted and propagated onto the bordering face /. Let 
and e-i be the other two edges of / such that v = €2 n e^. Propagating the sibling rays r and r' = S(£ (r, e±), v), 
where r' G Tr(v), involves computing the direction of these rays from face /' to /, which is accomplished using 



Snell's refraction formula. Suppose one of the sibling rays strikes e2 and the other strikes e^. Then the set of 
rays R(e\,w) is split into two sets of rays R(e2, w) and R{e^,w) via a binary search procedure. Sibling rays in 
these two sets, say r\ and T2 are computed and an approximation of the distance from w to v is computed. Then 
the vertex w originates a new set of rays. We shall show that if the rays are sufficiently close enough, in angular 
distance, then the approximation to the shortest distance from w to v is bounded by a factor of 1 + e. 

Further, the rays, R(e±,w) may critically reflect at e\, in which case the critical point of entry, say p\, needs 
to be determined. This generates a critical edge segment, CS = (pi,v') where v' is an endpoint of e±. All rays 
in R(e\,w) that are incident to points on the line segment CS reflect back onto the face /. The critical point of 
entry and the critical edge segment each become the origin of the following two sets of rays. See Fig. |2] We let 
r\ be the ray in R{e\, w) that strikes e\ at p\ and let r be the ray adjacent to r\ incident to e\ and at less than the 
critical angle. 

1. Set 1: Steiner Rays. Let 6 be the angle which the refracted image of r, termed tr, makes with the normal 
at edge e in face /'. The angular space w.r.t. edge e in the range [0, § — 0] needs to be populated with rays 
to ensure an e approximation. To achieve this, Steiner rays are generated from p\ onto face /' such that 

(a) There is a steiner ray, r s \ originating at pi onto face /' that is parallel to tr. 

(b) There is a set of rays originating from p\ onto face /', the rays uniformly subdividing the angular space 
between r s \ (which originates at p\) and (pi, v'). 

2. Set 2: Critically Reflected Rays. The other set of rays, which are critically reflected, are generated using 
the points on the critical edge segment CS. These rays have the property that they are parallel and leave 
the edge e and enter back onto face / at an exit angle which is the critical angle 9 C (Lemma [2. II) . Thus 
the critical edge segment originates a set of parallel rays each separated by a weighted distance of 8, the 
siblings denning these sets of rays being easily computed. See Fig. [3] Note that, fortunately, since two 
critical reflections cannot occur simultaneously on a shortest path (Lemma [2.21) . a critically reflected ray 
will not generate any further critical points of entry. In addition to the reflected rays, the entire critical 
segment from the critical source, pi, to the endpoint v' is also included as a potential shortest path to v' via 
the critical source. 

The above description of Set 1 gives the following property. 

Property 2.1 Two adjacent rays in Tr(v) are either parallel or originate from the same node w € Tr(v) 



Bounding the number of traced rays 

The shortest path from the source s to the destination can be split into sub-paths where each sub-path can be 
classified into types of sup-paths: 

1. A Type-1 path starts at a vertex in V and only refracts along the path. 

2. A Type-2 path uses critical segments, where no two critical segments occur consecutively, i.e every section 
of the path starts at a vertex, refracts and possibly critically reflects (just once) before either reaching the 
destination or another vertex. 

As evident from Lemma I2l2l these two paths suffice. 

Our correctness proof is split into two parts. The first part shows how to closely approximate the Type-1 
paths, which do not critically reflect in-between. The second part shows how to approximate Type-2 paths, 
which critically reflect in-between. 

In order to establish that our algorithm approximates Type- 1 paths, we first establish a bound on the number 
of rays to be propagated from the source. 

Lemma 2.5 Let the angle between two successive rays in Tr(v), Vu € V U C, be less than or equal to f^r (^) n ■ 
Then a Type-1 shortest path from a vertex uto a vertex v can be e-approximated using sibling rays in 7r(v). Here 
K is a large constant, e 1 = n ^+LW) ' n * s ^ e num ^ er °f vertices, w is the minimum and W is the maximum 
weight over all faces, and L is the maximum length of any edge in V. 

Proof: Given a Type-1 path P from, a vertex, say u, to another vertex, say v, we show that there exists a Type-1 
path P' using rays in our discretization so that P' closely approximates P. This path can use Steiner rays also. 
We show that the rays originated from u and successive critical points of entry, considered by the algorithm, 
are sufficient to discretize the space and ensure an e-approximation to any Type-1 path from u. We do so by 



(a) establishing a bound on the divergence of adjacent rays as they propagate (b) establishing an error bound on 
using steiner rays to approximate rays in Tr{v) and (c) bounding the distance of a vertex from a set of rays in 
Tr(w), w a vertex or critical source and establishing errors in using a sequence of vertices Wi,i = 1 ... k and 
their corresponding Tr(wi) to determine the shortest path. 

(a) Divergence of rays: We first consider the section of the path with no intermediate vertices. We show a 
bound on the angle between adjacent rays so as to limit the divergence of the rays in Tr(u). Let two rays, 
say r',r", originating from a vertex u be separated by an angle 7, and traverse across regions with weights 
ot\, ct2, ■ ■ ■ , dp respectively. This would be the regions through which the section of the shortest path, under con- 
sideration, traverses. Let the edges intersected by the rays be e±, e%, ■ ■ ■ , e p -\ respectively. W.l.o.g. we consider 
the case when the angle between the rays r' and r", is increasing due to refractions at these edges and for this 
case to occur we assume that W = ai > ai > . . . > a p = w > 1. Let 9[, 9' 2 , . . . , 9'_i be the angles at which 
the ray r' incidents on edges ex, €2, . . . , e p _i respectively. Let 62,63, ... ,9 P be the angles at which the ray r' 
refracts at edges e±,e2, ■ ■ ■ , e p _i respectively. Similarly, let 6[ + Si, 6' 2 + 82, ■ ■ ■ , x + Sp-i be the angles at 
which the ray r" incidents on edges e\, 62, ■ ■ ■ , e p _i respectively. And, let 62 + 62,63 + 63, . . . ,6 P + 8 P be the 
angles at which the ray r" refracts at edges e±, e2, ■ ■ ■ , e p _i respectively. Assume that 6 V < 6 cr iu ca i < ^; and, 
Vj, 5i is small. 

For every integer i € [l,p — 1], 
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For every integer i € [l,p — 1], 

ai sin (#■ + Si) = a i+ i sin (6 i+ i + 5 i+ i) 

aj sin 6 i cos Si + cos i sin Si = Qj + i sin 6i + i cos Si + i + a^+i cos 6i + \ sin Si + i 

Since S is small (—0.05 < S < 0.05), approximating sin<5 with S and cos S with one, the above equation 
becomes 

oti sin 6'i + aiSi cos 6\ ps a i+ i sin 6 i+ i + a i+ iS i+ i cos 6 i+ i 
aiSi cos 6^ « a i+ iS i+ i cos 6 i+1 (from^)) 
ai cos 6', „ 



■>i+i 



a i+ i cos V i+ i 



(2) 



Therefore, 
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ttt^-) = 0. Then the above becomes, 
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s P < —n 
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COS 0' 



From ([T]), aj > a^i, and 0- < 6i + \ we know that cos^ > cos6i + \. Since > 1, c 
Si+i > Si. Hence, S p > S p ^i > ... > Si. 



(3) 

> 1; from ©, 



We next consider the distance between two points such that the points are (i) at equal distance from a vertex 
v (ii) lie on adjacent rays r' and r" in Tr(v). Consider an edge on which the ray r' incidents at pi with an 
angle of incidence and refracts from pi with an angle of refraction 6i + i whereas the ray r" incidents at p' 9 
with an angle of incidence 6\ + Si and refracts from p' 2 with an angle of refraction 6i + \ + Si + i. See Fig. g] 
Consider two points pi,p2 located on rays r',r" respectively such that pi,p2 are in the region rj, with weight 




Note: All distances in the above figure are unweighted 



Figure 4: Refraction of two successive rays from a vertex 

cti, the weighted distance from u to either of these points is d, and HP1P2II — d\. Let r" be incident to e, at a 
larger angle than r'. Also, consider two points p%,P4 located on rays r', r" respectively such that Pz,Pa are in the 
region r^+i with weight cti+i, the weighted distance from u to either of these points being d', and ||j?3P4|| = d%. 
We wish to establish a bound on d2- Note that it suffices to restrict attention to the above case when both p^ and 
j>4 are in r^+i. If one of the points, say p4, is in ri then the bound on c?2 will be smaller since aii+i < cti. Also, 
let HP2P2II = \\P2PiW = d-2- We choose d" such that d'{ sin (0i+i + <5j+i) = sin (sin 6^ + Si). 

(d - d) = a i d 1 + a i+1 d 2 = . 7^—T\ a *+i«2 = ^ «i+id 2 ( 4 ) 

sin + dj J a i+ i 

We assume that all points at a weighted distance of d satisfy 

ccjdi < de' (5) 

for e' < e and determine conditions such that 

ai+1^2 < d'e' (6) 

a i+ id 2 

< aj + i((xi + X2) + y) (from triangle inequality) 

= sin7-^^sin0j + i+d / 1 sin (0< + 5j)+d 2 sin (0 i+ i + ^+i)| + |^^ cos 6> i+ i-d 2 cos (0 i+ i +^ + i)|) 

= |ai+idi sin7 - (d' - d) sin^+i + a{ + id[ sin (6>- + Si) + a i+ id 2 sin (0 i+1 + <$i+i)| + |(d' — d) cos - 
aj+id 2 cos (9 i+ i + 

< sin 7 — {d! — d) sin 9i + \ + (ajd^ + aj + id 2 ) max{sin {9[ + <5j), sin + <5j+i)}| + \ (d! — d) cos 0j+i — 
a i+ id 2 cos (0 i+ i + (for a* > a i+ i) 

= jajdi sin 7- (d' - d) sin# i+ i + (d' - d) sin + <5i+i)}| + |(d' - d) cos^+i -a i+ id 2 cos + <5i+l)| 

< jajdi + {d! - d)(sin(0 i+ i + <5 i+ i) - sin0»+i)| + |(d' - d)(cos#j+i - cos + <S»+i))| (from ©, ©, and, 
for < {ffi + Si) < (8 l+1 + 5 i+1 ) < tt/2, sin (Q[ + Si) < sin (6 i+1 + S i+1 )) 

< \de' + (d' — d)5i+i \ + \ — (df — d)Si+\\ (from ©, and, for small 
<de' + (d'-d)(2^ +1 ) 



Figure 5: Bounding the error when an optimal path traverses between two successive rays 



Since 6i+\ < S p , and, from ©, 2<5j+i < 25 p < 2^ L /3 p 8i < 2^ L 8 n2 5i (in the worst-case, p the number of 
edges intersected by a ray is 0(n 2 )) 



To satisfy ©, we need to have, 2^(3 n2 5i < e' i.e., ^ < * ) n V 



We thus get a bound on 7, 



7 = *<£(i)"V 



to achieve abound on the divergence of adjacent rays. 



(b)Error in using steiner rays: The above analysis will be used to bound the approximation of the distance 
from v to a point that lies between the two adjacent rays. However these two adjacent rays may not exist since 
a ray striking at the critical angle reflects back into the face. For faces fa, fi+i, let the edge ej = fi fl fi+i- For 

a ray which incidents on e.j at 9[ from f. L and refracts at 9{ along /i+i, let 3 = c ^°g . Suppose the ray r.; is 

incident on edge e, at point p with an angle of incidence greater than or equal to c j. Then the algorithm initiates 
a set 5 of Steiner rays from p. Among all rays in S, suppose the ray r 6 S makes the largest angle 9 r i with the 
normal to e, . Then B = — c ?l d % < . \ < 4-, for a large constant K and e 1 as defined below. 

Consider two successive rays r.; and rj from the same point, say v\. See Fig. [5] Let the ray r{ be incident 
at point p'l located on edge e, at an angle less than the critical angle. And, let the ray rj be incident at point 
pi located on edge e, at an angle greater than or equal to the critical angle. Let dp(x, y) represent the distance 
between points x and y along path P. Suppose an optimal shortest path OPT intersects edge e between p'[ and 
p\ at p[. Let OPT refracts at point p' x at an angle 9'. Further, suppose the ray r € S makes an angle 9" with 
the positive y-axis and suppose that for any ray r' € 5, the absolute difference in the angle r' makes with the 
positive y-axis and 9' is greater than or equal to \9' — 9"\. Consider a path P along the ray rj from v\ to p\ and 
along the ray r. Then the recurrence dp{v\,V2) = dp(vi,p\) + dp{p\,V2) represents the distance from point v\ 
to point V2 along path P. Let p'[ be at weighted distance d'[ from v± and, let p\ be at weighted distance d\ from 
v\. For d' = max(di,d' 1 / ) and a e being the weight of edge e and with a-id\ < de', we have a e ||piPill < d'e' 
(by the choice of 7 in equation ©). Then dp(v\,V2) < (1 + e')(dopr(vi,p' 1 ) + dp(pi,V2))- Let d be the 
optimal shortest distance from v\ to V2- Since there can be at most 0(n 2 ) critical points of entry along any path, 

expanding the recurrence yields dp{v\,V2) < (l + e') Cl?1 d. Let the error in computing the shortest distance from 
v\ to V2 be denoted by y. Then for small e', 

3^ < d Cl en 2 (8) 



(c) Error in reaching the destination: . We finally show that the distance from an origin u to another vertex v 
can be approximated by the sibling rays in Tr{u). We let £ be the sequence of edges traversed by the optimum 
path from u to v. Let the vertex v = e, n ej. Consider the sibling rays in the sequence of rays that intersect the 
sequence of edges £. Let the last edge in the sequence be e& before the shortest path reaches v after which the 
shortest path enters the triangle e,, ej, where v is common to ej and e&. We will compare the distance between 
the ray representing the optimal path from u to v and the rays in Tr{u). More specifically we will consider the 



Figure 6: Bounding the error when a vertex falls in between two successive rays 



rays in Tr{u) that intersect the sequence 8. We consider the two sibling rays n € Tr{u) and rj € Tr{u) such 
that the ray r,- L is incident to edge at pi, closest to v and the ray rj is incident to edge ej at pj, again closest to 
v. At least one of these sibling rays exist from the set of rays that intersect £. Note that the rays in Tr(u) that 
intersect the sequence £ split and intersect either both ej and ej or at least one of the edges ej or e^. See Fig. [6] 
We consider the first case. The second case is similar. Suppose \\vpi\\ = k, \\vpj\\ = lj. Letpi be at weighted 
distance di from u and pj be at weighted distance dj from u. For a = max(e?j, dj), from the bound (O we know 
that HpiPjH < d'e'. Let p\ be located on ray r; extended such that = k- Similarly, letpj be located on 

ray rj extended such that = lj. From above, we know that < max((dj + a ei li)e', (dj + a ej lj)e') 

where a ei is the weight along the edge ej and a ej is the weight along the edge ej. Let L be the maximum edge 
length and, let W be the maximum weight. Then ||f/jPj-|| < max^e', dje') + LWe'. Then at the vertex v, the 
error, X , i.e. the distance from the vertex to a point on a traced ray, is upper bounded as follows: 

X < LWe (9) 

Even in the case when the set of rays do not split and there is only one sibling ray closest to v to consider, the 
vertex v on the ray corresponding to the shortest path from u to v is at distance less than d w (u,v)e' from the 
sibling ray (from [7]). The rest of the proof is similar to the first case. 

From ([8]> and ©, and, given that there are C2n vertices on a shortest path where successive errors can accu- 
lumuate, the total error in computing the distance from u to v is at most dcne'(n 2 + LW), for c = max(ci, C2). 
Since this should not exceed de, we choose e' = n ( n -2+ LW ^ ■ Note that the approximation is achieved by rays 

which are closest to the vertex v. These rays are sibling rays in Tr{u). □ 



Lemma 2.6 Let P be a Type-2 shortest path from a vertex v to another vertex w onV with at least one critical 
segment CS in-between. If the angle between successive rays is as specified in Lemma |231 then an e-approximate 
Type-2 shortest path can be found from sibling rays in lZ(v),\/v G V |J C. 

Proof: The optimum path P can be partitioned into sub-paths, each sub-path going from a vertex v to a vertex w. 
Let this sub-path use a critical segment and be partitioned as: a path Pi from v to the critical segment, CS € /C, 
and a path P2 from CS to w. Let e be the edge on a face / such that the critical segment CS lies on e and reflects 
rays back onto face /. Also, let the critical segment CS have a critical point of entry cs. We first show that a 
good approximation to the path from CS to w can be found. If w is the endpoint of CS then we are done since the 
distance to the endpoint from the critical source is included in consideration. Otherwise consider rays generated 
from CS that are parallel and separated by a weighted distance of 5. As in the previous lemma, we consider in 
detail the case when w lies in-between two adjacent (parallel) rays from CS, T\ and r-i, which are sibling rays. 
As in lemma [231 the error in computing the distance from wtoa point on one of the rays is bounded by LWe 1 . 



The optimal path P2 has a first segment that is parallel to and lies in between a pair of rays r' and r". The source 
of the optimum path lies in between the origin of the two rays, origin(r') and origin(r"), that lie on CS. Since 
the weighted distance between origin(r') and origin(r") is at most 5, the error in determining a path along one 
of these rays is 5 + LWe 1 . Further the error in determing the shortest path from v to cs is at most y < dcie'n 2 
(equation. [8]). Combining the above errors, choosing 5 < e' and noting that there are at most 0(n) vertices on 
the shortest path gives the required error bound as in lemma [231 

□ 



3 Algorithm Details 

In the algorithm, the main procedure is the event handler. Event handling is categorized into the following types: 

1 . Shooting rays from a vertex 

2. Refracting a ray 

3. Finding the critical point of entry 

4. Shooting critically reflected rays from a critical segment and shooting Steiner rays from a critical point of 
entry 

The description of procedures that process the events is described below. The determination and handling of 
these events is further detailed with the pseudo-code. 



3.1 Main Procedure 

The main procedure is weightedEuclideanSP. It implements a version of continuous Dijkstra's algorithm as 
follows: An event heap is maintained, the heap comprising points ordered by their distance to the source vertex. 
The event may be a vertex, a point of strike of a ray, r, on edge or a critical source. Depending on whether the 
ray strikes the edge at an angle less than or greater than the critical angle, the ray is refracted or dropped from 
further consideration respectively. Let a ray and its sibling be such that one of them strikes e-i at an angle greater 
than the critical angle and the other at an angle less than the critical angle. The point at which a ray, starting from 
either of the sources of these rays, strikes the edge at the critical angle is determined. The distance to this point 
is determined and added to the heap. This point is termed as the critical source. Let e be an edge whose one 
endpoint is t. The algorithm terminates whenever the event is a vertex event with the vertex being t. 
The actions required at each event are listed below: 

(i) A vertex event triggers the procedure initiateVertexSource which initiates rays from a given vertex v. 

(ii) An event which represents a critical point of entry p\ uses the procedure initiateCriticalSource to initiate 
both critically reflected rays from the corresponding critical segment and Steiner rays from p\ respectively. 
From each critical segment, initially only two critically reflected rays are progressed as siblings. 

(iii) An event where a ray strikes an edge is handled in two parts. The procedure refractRay extends the ray 
across the current face /. It determines whether the given ray r' and its current sibling are incident onto the 
same edge of /. If not, it finds the correct siblings for both of these rays. Then it adds events corresponding 
to the rays striking the edges of the face /. The procedure findCriticaWointOfEntry finds the ray which 
originates from a vertex/critical point of entry and is incident to an edge at a critical angle. 

To improve the efficiency, as outlined before, we only choose to maintain sibling rays, rather than tracing all 
the rays that have originated from a point. Since an intermediate ray r between two sibling rays always has the 
same intersected edge sequence £(r) as one of these sibling rays, the intermediate ray is not required to be traced 
until it is necessary. These binary searches are detailed in the following procedures. 



Algorithm 1 weightedEuclideanSP 



1: Initialize distance between all pairs of vertices, d w (v, w), v, w G V to be oo. 

2: Push source s to the event queue 

3: while the event queue is not empty do 

4: evt <— pop the event queue 

5: {The following cases are both mutually exclusive and exhaustive} 
6: if evt represents a vertex v ^ t then 
7: initiate VertexSource(u) 
8: end if 

9: if evt represents a critical point of entry cs then 
10: initiateCriticalSource(cs) 
ll: end if 

12: if evt represents a ray r' which incidents to edge e, = fi n /j+i from face /, at an angle less than 
arcsinf-^-) then 

u 

13: refractRay(r', ej) 
14: end if 

15: if evt represents a ray r' with r 1 .sibling = r" such that r' incidents to edge = fi n /i+i from face /j 
at an angle < c j and r" incidents on the edge e, from face /j at an angle > 6 c i, for 6 c i = arcsin( ) 

fi 

(note that here r' and r" can be replaced one for the other) then 
16: if source of ray is not on a critical segment then 
17: nndCriticaiPointofEntryCr' , r" , e{) 

18: end if 
19: end if 

20: if evt corresponds to a vertex event where the vertex is t then 
21: empty event queue 
22: end if 
23: end while 



3.2 Procedure initiate VertexSource 



This procedure initiates rays from a vertex v, including the source. Every face that includes the vertex v is 
considered. For each such face /, and the edge ej = (u, w) G / that is not incident to v, two rays from R(e{,v), 
r • and r'( are determined such that rj is closest to u and r • closest to w (in fact the procedure discovers successive 
rays in Tr(v) that are closest to the endpoints of the edge ej). These rays and r- are the sibling rays which will 
be propagated further. The distance at which these rays will intersect ej is determined and added to the heap. 

A binary search on Tr{v) is used to ensure that rays in Tr(v) are correctly generated along with sibling rays. 
The following is clear from the above discussion: 



Algorithm 2 initiate VertexSource(t>) 

Require: A vertex v. Let F be the set comprising all the faces incident to v. Let the set of edges defining faces 

in F be E' . Also, let E"(C. E') be the set of edges which incident to v. See Fig. [TJ 
Ensure: Initiates rays from the vertex v 

1: for Each vertex v ' ^ v such that v' incident to a face / € F do 

2: Binary search over the rays originating from v to determine two successive rays, 7$, r'', such that 
for ej,ej EE' — E",v' = ej n ej, where r\ = arg mm. r&R ^ ei ^{d(p(r), v')} and r'- = 
argmin re ^( e „){ci(p(r), v')}. The rays r[, r'j are the ones traced further. 

3: end for 

4: for Each edge ej in E' — E" do 

5: Let rays r-, r" be sibling rays in R(ei,v). Let the point p\ = r[ n e, and, let = r-' n ej. Also, let / 
be the face on which the line segments vp[ and vp" resides. Push new events to the heap representing the 
intersection of ray r\ with ej, which occurs at a distance dy = a/||«j^|| from v and the intersection of ray 
r'l with ej, which occurs at a distance dj« = a/||vpj'|| from u. 

6: Let f ' and v" be the endpoints of ej closest to p\ and p", respectively. Update distance from source v to v' 
in heap by min{d w (v, v'), dy + d(p' { , v')} and to v" by min{(i TO (t>, w"), dj» + d(pi>',v")}. 

7: end for 



Lemma 3.1 Procedure initiateVertexSource correctly initiates the tree Tr(v). 
3.3 Procedure initiateCriticalSource 

When a critical point p on edge e = / n /' of face / is encountered, there are two classes of rays that need to 
be propagated. One is the set of rays that are to be critically reflected back onto /. And the other is the set of 
Steiner rays over the face /' that are required to be generated from the critical source p. As required in the proof 
of Lemma 1231 these Steiner rays are required to achieve the desired approximation as they serve to fill in for the 
set of rays on face /' that strike the edge e at an angle just less than the critical angle. More specifically, let n, 
with r2 as the designated sibling, be the ray that strikes e at point p. Let this ray r\ be an element of Tr(w) and 
let u be the origin of r\. Note that u may be w itself or a critical source. Consider the ray r% which is adjacent to 
r\ in Tr(w) and which would be refracted at edge e. Simply refracting r% would not be sufficient since the rays 
incident to e at angles between the angle a at which r% strikes e and the critical angle, refract to cover a wide 
range of angles at the adjacent face. These rays are replaced by rays that originate from p. A node corresponding 
to the critical source p is created as a child of u in Tr(w) and a set of rays TZ(p) are initiated. Sibling rays that 
are propagated forward are computed as follows: One sibling ray, termed r\, is the ray almost parallel to e, in 
fact which makes an angle e'/K w.r.t. e and the other is the sibling of the ray n, termed t2, which struck e. 

These rays are propagated across the face /'. Let e2 and e\ be the edges of /' that r% and strike. If e2 
is the same as then the rays r2 and are siblings. Otherwise, new siblings are created using the procedure 
findSplitRays. Let v be the vertex common to e2 and e^. The shortest distance to v from amongst the rays in 
Tr(w) lying between T2 and r$ is computed and the event heap corresponding to the distance from the source 
to the vertex source v is updated. Note that when the shortest distance to v will be discovered as the minimum 
distance event in the heap, v will be a vertex source. 

Lemma 3.2 Procedure initiateCriticalSource correctly computes a new node w corresponding to a critical 
source on an edge e and the set of rays lZ(w). It also correctly computes sibling rays incident to edges other than 
eoff. 



Algorithm 3 initiateCriticalSource(pi ) 

Require: Critical point of entry p\ on edge e = (u, v) and t\ € Tr{w) the ray which caused the critical point of 
entry p\. Critical segment seg = (pi,v). Let / be the face onto which the critically reflected rays from seg 
traverse. Let 9 C be the critical angle of incidence of edge e. Also, let T\. sibling = 7"2, and, let P2 = ?~2 fl e. 
See Fig. E 

Ensure: Initiates critically reflected and Steiner rays from seg and p\ respectively 

1: Let pi be a point at \\p\v \\ — 5' distance from p\, for 5' < e'. Let r', r" be critically reflected (parallel) rays, 
making an angle — 9 C with the normal to e, and originate from points p\ and p' x respectively. Let p' be the 
point at which the ray r' incidents on an edge e' of / and, let p" be the point at which the ray r" incidents 
on an edge e" of /. Note that e! is not necessarily distinct from e". Set r', r" as siblings. Push new events 
evt', evt" to the heap which represent the intersection of r', r" with the edges of / where evt' occurs at a 
distance a/Hpip'H frompi and evt" occurs at a distance a e ||pip'i || + "/IIp'iP"!! from p\. Also, push an event 
to the heap which corresponds to the critically reflected ray reaching vertex v. This event occurs at a distance 
a e ||pi«|| frompi. 

2: if e' / e" then 

3: findCriticallyReflectedray(r' , r" , v') where v' is the vertex common to e' and e" 

4: Let r3 be the ray adjacent to v\ in Tr{w) which refracts when it incidents via / to a point P3 located on 
line segment p\P2- Suppose 7-3 refracts from p% at an angle 6 onto face /'. A set S of Steiner rays, which 
make angles in the range [j^, ^ — 9] with edge e and traverse over /', are added to Tr(w), at a new node 
pi, i.e. the rays originate from p\. Here, K is a large constant. Note that if | — 9 < then the range 
[jr , f — 9] is 0, and no rays originate from p\. Let S S be the ray which makes an angle with edge 
e. Suppose the ray r± first intersects with an edge e^ of /' at point p,\. Set the rays 7-2 and r<± as siblings. 
Push new event to the heap which represents the ray 7-4 traversing to p^ and this event occurs at distance 
a f \\PiP<i\\ from p\. Trace T2 over face /' by pushing an event corresponding to the intersection of T2 with 
an edge 62(7^ e) of /', termed p' 2 . Note that e2 is not necessaiily distinct from e±. 

5: if e2 = e4 then 

6: Update shortest distance to endpoints of e2 from p' 2 and p±. 
7: else 

8: findSplitRaysO and update shortest distance to the endpoints of e2 and e^ from the siblings. 
9: end if 
10: end if 



3.4 Procedure refractRay 

In procedure refractRay, a subset of rays in Tr{s) that strike an edge a before traversing a face / are refracted. 
These rays lie in between r' and its sibling r" that strike ej. Let ej, e 3 ■, e^ be the edges of /. And let v be common 
point to ej and e^. 

Suppose r and r' , when refracted, strike different edges of the face, say ej and e^, respectively. If r and 
r' have originated from a critical segment source, CS then a ray r c originating from CS and that strikes v is 
computed. Also are determined the siblings of r and r', which would be the rays adjacent to r c in Tr{CS). 

Alternatively, r and r' originated from a vertex source w. Then the ray in Tr(w) that strikes v is computed 
and the siblings of r' and r", n and r% respectively, are determined. This is done using procedure findSplit 
Rays. The procedure findSplitRays proceeds as follows: The critical ancestor paths, P' and P" of r' and r" are 
determined. The siblings must lie in the set of rays generated from vertices along the paths P' and P". These set 
of rays are searched via binary search. The events at which these rays strike a and ej, respectively, are computed. 

Lemma 3.3 Procedure refractRay correctly processes refracted rays across a face f and updates sibling rays 
incident onto edges of the face f. 



Algorithm 4 refractRay(r' , ej) 

Require: The ray r' which incidents on edge ej at an angle less than 6> c j. Let r'. sibling = r" on the same edge 

Ensure: The rays r' and r'. sibling are traced further. 

1: Let / comprising edges ej, ej, e^ be the face over which the rays r' and r" are to be traced 

2: if r' and r" traced forward are incident onto different edges ej and e^ then 

3: if r' and r" are originated from the same critical segment then 

4: fmdCriticallyReflectedRayfr' , r" , v ) where v = ej n e&. 

5: else 

6: findSplitRays(r / , r", ej, e^) 

7: end if 

8: else 

9: Push events corresponding to rays r' and r" striking edges ej and e^. 

10: end if 



Algorithm 5 findSplitRays(r' , r", ej, ej) 

Require: For the rays r' and r'. sibling = r", such that r' , r" € Tr{w), let v' be the least common ancestor of 
r' and r" . Let P' be the critical ancestor path of ray r' with v' as the origin. Let r' be incident to edge e^ at 
point p' and r" be incident to edge ej at point p", for v = e^ C\ ej. Similarly, let P" be the critical ancestor 
path of ray r" with v' as the origin. Let REG be the region bounded by the rays in the critical ancestor path 
of r' and the critical ancestor path of r", rays r' , r", and the segment p'p". Let Ri be the set of rays such that 
a ray r G R\ if and only if r originates from a critical point of entry/vertex located on the critical ancestor 
path P' and ray r lies in the open region REG. Similarly let R2 be the set of rays such that a ray r € R2 if 
and only if r originates from a critical point of entry/vertex located on the critical ancestor path P" and ray 
r lies in the open region REG. See Fig. Ha). 

Ensure: Finds two rays n, r2 which intersect edges ej, ej respectively such that among all the rays in R± and 
i?2» the rays n, r2 are the closest rays to vertex u along the edges ej and ej respectively. 

1: Do binary search of the rays in Ri and the rays in R2 to find two rays n and T2 where n intersects ej and 
r2 intersects ej such that n and T2 are either successive rays originating from the same origin or adjacent 
origins on a critical ancestor path. Set r\,r' and r", T2 as siblings. Push two events to the heap where one 
represents the intersection of ray n with ej and the other represents the interaction of ray T2 with ej. 

2: Let rays n, r% intersect edges ej, ej respectively at points p\,p2- Thenmm(d w (s,pi)+a e d(piv),d w (s,p2)+ 
a e d(p2, v), is the event required to be update the shortest path to v in the event heap. 



3.5 Procedure findCriticallyReflectedRay 



This procedure is invoked in the previous procedure when the ray to be refracted arises from a critical segment 
CS. Consider when two critically reflected sibling rays r', r" from the same critical segment CS are propagated 
across a face and intersect two distinct edges ej , ej respectively. This procedure finds new rays r' sib and r" sib which 
are siblings to r', r" respectively so that r' sib and r" ib originate from CS. 



Algorithm 6 fmdCriticallyReflectedRay(r / , r" , v) 

Require: Two critically reflected rays r',r" originating from the same critical segment, say CS, with 
r'. sibling = r"; r' incidents on edge e, and r" incidents on edge ej for v = ej H ej. See Fig. [3] 

Ensure: Finds a critically reflected ray r originated from a point on CS such that r incidents to v. Also, finds 
new siblings r' sib and r" ib for r' and r". 

1: Let p' = origin(r') and p" = origin(r"). Let / be the length of the edge sequence intersected by r' i.e., 
/ = Consider a critically reflected ray r such that r = origin(p) on the line segment p'p" which 

passes through v, is parallel to r' and £(r) is the same as the first / — 1 elements of £(r r ). Update event in 
heap corresponding to d w (CS, v), the weighted distance from CS to v. 

2: Construct the ray r' sib , r' sib = origin(p' sib ) where p' sib is located on the line segment p'p" such that 

(i) r' sib is parallel to r' 

(ii) £{r') is the same as £(r' sib ) 

(iii) d{r,r' sib ) = 6 (for S < e'). 

And construct the ray r" sib be the mirror image of r' sib w.r.t. r, i.e. r" sib = origin{p" sib ) where p" sib is located 
on the line segment p'p" such that 

(i) r" ib is parallel to r" 

(ii) £{r") is the same as £{r" ib ) 

(iii) d(r, r" ib ) = 8 (for S < e'). 

Set r', r' sib and r", r" sib as siblings and compute current estimates of the shortest path to the endpoints of 
and ej, via the sibling rays. 

3: Push events to the heap which represents the weighted distance at which ray r' sib intersects ej and the ray r" sib 
intersects ej. 



3.6 Procedure findCriticalPointOfEntry 

While propagating a ray it may happen that r' and r" are two sibling rays that originate from the same point and r' 
is incident to ej at an angle less than the critical angle while r" incidents to ej at an angle greater than the critical 
angle. Procedure findCriticalPointOfEntry handles this case and uses the procedures wrapperToFindHitAtAngle 
and Find-Hit-At-Angle together to find the ray r which incidents at the critical angle on edge ej such that the ray 
r traverses across the same edge sequence as r' and r originates from v. 

Since a vertex must exist between two critical points along a shortest path, we do not progress a ray r striking 
an edge at an angle greater than or equal to critical angle whenever r is originated from a critical segment (see 
procedure wrapperToFindHitAtAngle). For the sake of completion, we have listed the procedure Find-Hit-At- 
Angle from along with wrapperToFindHitAtAngle'm the Appendix. 

4 Correctness and Analysis 

Lemma 4.1 The algorithm correctly determines sibling rays in Tr(u) for every source u, and computes an 
e- approximation to a Type-1 shortest path from u to v for any two sources u, v. 

Proof: Suppose that there exists a Type 1 shortest path between u and v , traversing a sequence of edges £. 
Lemma 1231 shows that this path can be found using sibling rays. We show that the algorithm maintains sibling 
rays for edges in £ that enables the determinations of traced rays such that these rays can be used to find an e- 
approximate shortest path from u to v. First consider the case when the shortest path lies in between two sequence 
£ . These adjacent rays will be found from sibling rays by the algorithm as shown below. The second case occurs 
when the shortest path lies outside the set of rays, R. In the second case, the vertex v must be less than d w {v)e' 
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comprise r and r 3 



(a) Sibling rays are separated by a vertex v 
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critical ancestor path of r' 
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(b) One sibling ray incidents at an angle less than the critical angle and the other incidents at an angle 
greater than or equal to the critical angle on a 



Figure 7: Finding sibling rays using critical ancestor paths 



Figure 8: Finding the ray which incidents at critical angle on edge 



Algorithm 7 findCriticalPointOfEntry(r / , ej) 



Require: Ray r' which incidents at point p' located on edge ej = fa D /i+i with an angle of incidence less than 
or equal to critical, and, ray r'. sibling = r" which incidents at point p" of edge with an angle of incidence 
greater than or equal to critical. See Fig. [7Jb). 

Ensure: Finds a ray that originates from a critical source/vertex and incidents to a point p located on the line 
segment p'p" at an angle equal to the critical angle. 

1: Let the sets R\ and R2 be denned as in procedure findSplitRays. Do binary search on the critical sources 
defining the rays in R\ and in R2 and subsequently on the rays to find two rays r\ and T2 where r\ incidents 
on ej at an angle less than or equal to critical and T2 incidents on a at an angle greater than or equal to 
critical. 

2: if original) = origin^) then 

3: wrapperToFindHitAtAngle(ri , r-i, ej, origin{r\)) 

4: else 

5: Suppose origin{r2) is an intermediate vertex on the critical ancestor path of r" (the other case is symmet- 
ric). And, let r be the ray with origin same as origin{r\) that led to the critical source, originfo). Then 
we invoke wrapperToFindHitAtAngle(ri , r, e^, origin(r{)) 

6: end if 



from a ray (in fact a traced sibling ray) in R, since two points at distance d on adjacent rays, are separated by a 
distance of atmost de 1 (due to|7]in Lemma |23T ). The algorithm computes the estimate of the shortest distance to 
a vertex v whenever sibling rays are computed at every edge incident to v. 

It thus suffices to show that sibling rays are computed correctly. Initially the siblings are computed correctly in 
procedure initiateVertexSource (Lemma [3TTb and each ray is refracted correctly in procedure refractRay (Lemma 
13.31) . We consider the other procedures whose correctness is not immediate from their descriptions. Consider 
the procedure findSplitRays which computes the siblings r\ and T2 closest to a vertex v at some stage in the 
algorithm. The rays n and T2 are thus adjacent rays and the distance to v via the rays n and T2 is updated. 
Suppose both the rays r\ and T2 have the same origin. Recall that Procedure findSplitRays finds the siblings 
and due to ((TJ) and ©, an e-approximation to the shortest distance from origin(r') to v is ensured. Otherwise, 
origin(ri) ^ originfa)- In this case procedure findSplitRays finds the siblings via a binary search on the rays 
in the critical ancestor paths of two siblings, say r' and r" such that rays r\ and are part of the set of rays lying 
in between r' and r". Let r be the ray which initiated originfa). And let ray r has as its origin p = origin(n). 
Also, we know that the ray r\ is a successor to ray r since it is adjacent to r%. Let originir^) be located on 
an edge e'. Suppose n intersects ef at point p\. See Fig. Ha). As described , by construction (procedure 
initiateCriticalSource and Lemma [372b . ray ri is parallel to r2, starting at p\. Due to ©, the distance between 
pi and originfo) is bounded. These facts, together with Q, guarantee the e-approximation in computing the 
shortest distance from origin(r') to v. 

Note that initiateCriticalSource requires the procedure findCriticalPointOfEntry whose correctness is clear 
from the description. □ 

Lemma 4.2 Let P be a Type-2 shortest path from a vertex v to another vertex w on V with a critical segment 
CS in-between. Then the algorithm determines a pair of traced rays in 7Z(v) that can approximate the shortest 
path, P. 

Proof: The proof of lemma 1231 shows that sibling rays in Tr(v) and sibling rays originating from CS will 
determine the approximate shortest path. The optimum path P can be partitioned into two, a path Pi from v 
to CS £ JC and a path P2 from CS to w. Let e be the edge on a face / such that the critical segment CS lies 
on e and reflects rays back onto face /. Also, let the critical segment CS have a critical point of entry cs. We 
first show that a good approximation to the path from CS to w can be found by showing that sibling rays are 
correctly maintained. In fact, we need to find the sibling rays closest to vertex w when the rays from CS strikes 
edges incident to w. Note that initially, two rays from the end points of CS are generated onto face / correctly 
in procedure initiateCriticalSources (Lemma 13^21 . as the sibling rays. Subsequently, as these rays strike edges of 
the faces they are refracted and siblings updated.. If the rays strike different edges e± and e^, e\V\e2 = w, after 
striking an edge e' of a face comprising edges e', e\ and ei, the appropriate siblings are computed in procedure 
findCriticallyReflectedRay. Also the shortest distance path to w starting at the critical segment CS is computed. 
This path has a first segment that is parallel to the sibling rays, r' and r", that are initially generated in procedure 
initiateCriticalSource and maintained during the progress of the propagation. Furthermore it has its origin, 
p = origin(r) located in-between origin(r') and origin(r"). Such a path is found in Step 1 of procedure 
findCriticallyReflectedRay. If the sibling rays do not split but strike one edge, say e±, the distance via the closest 
sibling to w is computed as required in lemma 1231 

Finally, a Type-1 path between u and cs is found correctly by the algorithm as specified in lemma |4~T1 □ 

Theorem 4.1 The algorithm computes an e-approximation to an exact shortest path from stot in 0(n G lg(-p-) + 
n 4 lg(^)) time complexity. 

Proof: The correctness follows from the previous two lemmas and by splitting the optimal path from s to t into 
parts, each of which can be classified as either a Type 1 or a Type 2 path. 
We bound the time complexity by 

1. Computing the time required to determine sibling rays, i.e. rays that are traced by the algorithm in find- 
SplitRays. 

2. Computing the time required mfindCriticalPointofEntry. 

The other procedures called in weightedEuclideanSP will be shown to have lower order of complexity. The 
introduction of sibling rays is done in procedures findSplitRays and findCriticallyReflectedRay : When the algo- 
rithm determines new siblings, it is required to do binary search over lZ(v),v € C. The total number of rays in 



1Z(v),v € C are 0(^r(jr) n ). It takes 0(lg(;^(-p-) n )) steps to do binary search over these rays where each 
step involves tracing the ray forward and computing the intersection with the edge sequence to determine the 
point at which the ray strikes one of the edges ei or edge ej. Since an edge can appear in an edge sequence 0(n) 
times (Lemma |23l ), the length of an edge sequence is 0(n 2 ). This also implies that the size of the sets R± and 
i?2 (which are defined in procedure findSplitRays), are 0(n 2 ). Consider the critical ancestor path P such that the 
origin of P is v and P comprises the maximum possible number of critical points. Let S be the set consisting of 
all the critical points of entry in path P together with v. Note that if we do a binary search over the set of rays in 
lZ(v'),v' € S, we need not do binary search over 1Z(v"), v" € S. It then remains to determine the vertex v' £ S 
where binary search is required. This can be done by considering the vertices on the critical ancestor paths, which 
are of size 0(n 2 ). Therefore, the work involved for each invocation of findSplitRays is 0(n 2 lg(^p-(-p-) n )) and 

since this may be required for all pairs of vertices in V, the total work is 0(n 6 lg(-p-) + n 4 lg(^r)). 

With similar reasoning, the total work due to binary searches in the procedure findCriticalPointOfEntry takes 
0(n 6 lg(f) + n 4 lg(^))time. 

We now consider the time required for other procedures. The main algorithm also invokes the procedures 
initiateVertexSource, initiateCriticalSource and refractRay. The first two procedures use binary search to find 
sibling rays. This is done only O(l) times for every procedure call. And these procedure are called 0(n 3 ) times 
since there are 0(n 2 ) critical points of entry (Lemma |24l ) and there can be 0(n 3 ) possible combinations of vertex 
and critical points of entry. Next consider the time for procedure refractRay. Each invocation of this procedure 
results in processing a ray and its sibling. The complexity of the calls to procedure findSplitRays during the pro- 
cedure is detailed above. For the complexity of the remainder operations we note that the procedure refractRay 
requires O(l) work when called for each strike of sibling rays onto an edge of a face. The number of sibling rays 
is bounded by 0(n) for each source v since a pair of sibling rays can be associated with a vertex which caused 
the set of rays in 1Z(v) to split into two. This bounds the number of calls to refractRay to 0(n 4 ). Similarly, 
since there are 0(n 2 ) critical points of entry (Lemma [2~4l ). there can be 0(n 3 ) possible combinations of vertex 
and critical points of entry such that it is required to find the sibling rays in procedure findCriticallyReflectedRay. 
Hence the procedure Find-Hit-At-Angle is invoked 0(?i 3 ) times, whereas each invocation takes 0(n 2 ) time i.e., 
the total time spent in procedure Find-Hit-At-Angle during the entire algorithm is 0(?i 5 ). Since an edge sequence 
is of length 0(n 2 ), the work involved in tracing the ray incident to a vertex from a critical segment takes 0(n 2 ) 
time (line 2 of procedure findCriticallyReflectedRay). As there are 0{n) vertices and 0(n 2 ) critical segments 
(Lemma l24l ). the total work involved takes 0(n 5 ), which is subsumed by the complexity of other procedures. □ 



5 Conclusions 

This paper presents an algorithm with (9(n 6 lg(-p-) + n 4 lg(^)) time complexity to find the shortest path from 
s to t in a weighted subdivision. This algorithm is of better time complexity than the existing polynomial time 
algorithm and it is better than the pseudo-polynomial time algorithms whenever ^= > n 5 . 

References 

[1] P. Agarwal, S. Har-Peled, and M. Karia. Computing Approximate Shortest Paths on Convex Polytopes. 
Algorithmica, 33(2):227-242, 2002. 

[2] L. Aleksandrov, A. Maheshwari, and J.-R. Sack. Approximation Algorithms for Geometric Shortest Path 
Problems. In Proceedings of the ACM Symposium on Theory of Computing, pages 286-295, 2000. 

[3] L. Aleksandrov, A. Maheshwari, and J.-R. Sack. Determining Approximate Shortest Paths on Weighted 
Polyhedral Surfaces. Journal of the ACM, 52(l):25-53, 2005. 

[4] S.-W. Cheng, H.-S. Na, A. Vigneron, and Y. Wang. Approximate Shortest Paths in Anisotropic Regions. 
SIAM Journal on Computing, 38(3):802-824, 2008. 

[5] M. Lanthier, A. Maheshwari, and J.-R. Sack. Approximating Weighted Shortest Paths on Polyhedral Sur- 
faces. Algorithmica, 3 0(4): 5 27-5 62, 2001. 

[6] C. Mata and J. Mitchell. A New Algorithm for Computing Shortest Paths in Weighted Planar Subdivisions 
(Extended Abstract). In Proceedings of the ACM Symposium on Computational Geometry, pages 264—273, 
1997. 



[7] J.S.B. Mitchell and Christos H. Papadimitriou. The Weighted Region Problem: Finding Shortest Paths 
Through a Weighted Planar Subdivision. Journal of the ACM, 38(l):18-73, 1991. 

[8] Z. Sun and J. H. Reif. BUSHWHACK: An Approximation Algorithm for Minimal Paths Through Pseudo- 
Euclidean Spaces. In Proceedings of the International Symposium on Algorithms and Computation, pages 
160-171, 2001. 



[9] Z. Sun and J. H. Reif. On finding approximate optimal paths in weighted regions. Journal of Algorithms, 
58(1): 1-32, 2006. 



6 Appendix 



Algorithm 8 wrapperToFindHitAtAngle(r' , r", e^v) 

Require: Let the critical angle associated with edge e.j be 9 C . The ray r' incidents at point p' located on the edge 
&i = fi H /j+i with an angle less than or equal to 9 C and the ray r" incidents at point p" located on the edge 
ej with an angle greater than or equal to 9 C , for origin(r') = origin(r") = v. See Fig. [8] 

Ensure: This is a wrapper to the procedure Find-Hit- At- Angle. 

1: if n.sourcetype != critseg then 

2: Define the interval I as [p',p"]. 

3 : cs <- Find-Hit- At- Angle(# c , I , e , , v , r ' ) 

4: Let r c be the ray causing cs. Set n , r c as siblings. Let 9 be the angle at which r' refracts from p', and, let 
seg be the critical segment due to cs. Push new event corresponding to the critical point of entry cs into 
the heap, with seg, 9, 9 C as associated data. 

5: end if 



Algorithm 9 Find-Hi t- At- Angle(# c , I,ei,v,r') (procedure from Q) 

Require: an interval /, critical angle of reflection 9 C , edge such that ej n I = I, ray r' which intersects an 
endpoint of /, point v which is the origin of sibling rays incident to endpoints of I 

Ensure: Find the point in I (if there is one) such that the refraction path from the root of I through the last edge 
sequence will hit the point at the angle of incidence 9. If the function succeeds in finding the desired path, it 
returns the critical point of entry in the interval /; else, returns NIL. 

1: Starting with an angle of incidence 9 C at edge e$, determine the angles of incidence at each edge in the last 
edge sequence of r'. This can be done by following back pointers to predecessors of intervals and using 
the fact that the angle of incidence at e, uniquely specifies the angle of refraction from ej_i, which in turn 
uniquely specifies the angle of incidence at It may be impossible at some stage to achieve the desired 
angle of refraction, due to the constraint that the angle of incidence be less (in absolute value) than the critical 
angle; in this case, we stop and return NIL. This process will eventually give the angle j3 such that a refraction 
path that starts from v at an angle (3 will refract through the last edge sequence or r' and strike e« at an angle 
of incidence 9 C . 

2: Now do ray tracing from v at an angle /3, following the last edge sequence or r'. If the refraction path leaves 
the sequence, then stop and return NIL. Otherwise, the refraction path will eventually cross edge e, at some 
point x (and we know by construction that it will strike at point x with the desired angle of incidence 9 C ). If 
x G /, then return x; otherwise, return NIL. 



